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We study a Hamiltonian system describing a three- spin- 1/2 cluster-like interaction competing 
with an Ising-like anti-ferromagnetic interaction. We compute free energy, spin correlation func- 
tions and entanglement both in the ground and in thermal states. The model undergoes a quantum 
phase transition between an Ising phase with a nonvanishing magnetization and a cluster phase 
characterized by a string order. Any two- spin entanglement is found to vanish in both quantum 
phases because of a nontrivial correlation pattern. Neverthless, the residual multipartite entangle- 
ment is maximal in the cluster phase and dependent on the magnetization in the Ising phase. We 
study the block entropy at the critical point and calculate the central charge of the system, showing 
that the criticality of the system is beyond the Ising universality class. 

PACS numbers: 03.65.Ud, 05.30.Rt, 03.67.-a, 42.50.-p 



I. INTRODUCTION 

The interplay between quantum information and sta- 
tistical mechanics has given rise to a new trend in con- 
temporary physics research. On one hand, quantum in- 
formatics provides new views into statistical physics, with 
spin-offs that could lead to decisive progress in the field; 
on the other hand, strongly correlated quantum statisti- 
cal systems, considered as platforms for future quantum 
computers, naturally provide correlations that quantum 
information aims to exploit as a resource [U [2]. Cold- 
atom quantum simulators play an important role in this 
context because they provide quantum statistical sys- 
tems, beyond traditional condensed matter realizations 
[3]. Within this interdisciplinary field, quantitative anal- 
ysis of the entanglement encoded in a given state of the 
statistical system provides precious informations on its 
physical properties [IHS]. 

This article deals with a one dimensional statisti- 
cal system formulated in the cross-fertilization area de- 
scribed above. The reference system is provided by cold 
atoms in a triangular optical lattice [7]. For a suitable 
choice of the parameters, such a system can be consid- 
ered as a spin system with a specific ring-exchange inter- 
action in a triangular lattice that, in turn, can be recast 
into a "zig-zag chain". Remarkably, the ground state of 
the system is the so-called cluster state [8 . The setup 
provides a physical platform for the one-way route to 
quantum computation, where the algorithm consists in 
suitable measurements aiming at reconstructing the high 
degree of entanglement characterizing the cluster state 



[9]. Interestingly, besides the three-spin ring-exchange 
interaction, various two-spin interactions can be shown 
to emerge in the system. Therefore, the cluster interac- 
tion competes with the exchange one by tuning a control 
parameter [8j. The interplay between these two interac- 
tions has been independently considered in the context 
of quantum information, in order to estimate the effects 
of local perturbations on the cluster states [1QHT2]. It re- 
sults that the correlation pattern characterizing the clus- 
ter state is robust up to a critical value of the control 
parameter, meaningfully defining a "cluster phase"; for 
larger values of the control parameter the system is in 
an Ising phase. The two phases are separated by a conti- 
nous quantum phase transition (QPT) [13]. Interestingly 
enough, the cluster phase is characterized by a diverging 
range of localizable entanglement [14 that can be traced 
back to a non- vanishing so-called string order parameter 
115]. 

Similar findings were recently obtained for a slight vari- 
ation of the model discussed above, that we shall call 
cluster-Ising model (CIM), whose Hamiltonian reads [16] 

N N 

where af, a = x^y^z^ are the Pauli matrices and, ex- 
cept otherwise stated, we take periodic boundary con- 
ditions CF%_^]^ = cr^- We note that the nearest-neighbor 
and the next-to- nearest- neighbor exchange in the two- 
and three-spin interaction, respectively, involve orthogo- 
nal spin orientations (in contrast with the model studied 
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in [TQ1[TT]); therefore all the spin components participate 
to the interaction. 

In this article we intend to study the statistical me- 
chanical features of this system, by analyzing correlation 
functions, entanglement and quantum phase transitions. 
We shall focus in particular on the remarkable interplay 
between the competing "cluster" and "Ising" phases. 



II. OVERVIEW 

For A = the ground state is a cluster state de- 
fined as the unique common eigenvector of the set of 
commuting Hermitian operators (known as "stabilizers" ) 
Kj = crj_icr|crj+i: Kj \C) = +1 \C) [I?]. The quench of 
the cluster to antiferromagnetic order occurs through an 
exotic quantum phase transition from an antiferromag- 
net to a phase with a hidden order - the cluster state 
provides an example of Symmetry Protected Topological 
Order (SPTO) by a Z2 x Z2 symmetry [16, 18, 20^ (see 
[2T] for applications of the notion of SPTO to quantum 
computation). In the present article we shall corroborate 
these findings by studying the degeneracy of the ground 
state of the system through a mapping of the spins onto 
Majorana fermions (see Fig. |3|. For the cluster state 
A = the degeneracy of the ground state is due to four 
uncoupled Majorana fermions. 



We shall obtain in Sec. IV the partition function of the 
system together with the exact expressions of the correla- 
tion function i?J^(T) = {(jJaJ)T^ with 7 = x^y^z^ and T 
temperature, in terms of Toeplitz determinants. A non- 
trivial pattern of spin-spin correlation functions emerge: 
Rji{T)^ together with the magnetization along z identi- 
cally vanish; R^i{T) identically vanishes, unless \j — l\ is 
a multiple of 3; R^i{T) is nonvanishing. 

The cluster and the Ising phases are characterized by 
the string order parameter Oz and the staggered magne- 

Both will be ob- 
In particular, we 
1)8, showing that 



tization rriy, respectively. See F ig. [I 



tained exactly in Sees. |IV A| and [T 
find Oz ^ {I — A)^^^ and rriy ~ (A 
the quantum phase transition is not in the Ising univer- 
sality class (that would have yielded /3 = 1/8). 

The entanglement pattern will be investigated ir 
Sec. ~ 
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Two-spin entanglement identically vanishes, 
both in the "thermal" ground state and in the ground 
state with broken symmetry, both at T = and at any 
finite temperature. It is therefore unable to detect the 
quantum phase transition. Nevertheless residual (mul- 
tipartite) entanglement r is nonvanishing and is viewed 
as a global figure of merit. We find that it saturates 
at its maximal value r = 1 for any A in the "ther- 
mal" ground state. However, interestingly, when one 
looks at the symmetry-breaking ground state, multipar- 
tite entanglement saturates in the whole cluster phase, 
but it decreases for A > 1, marking the critical point with 

dxr{X) -cost X (A - 1)-!/^. See Fig.^ The mul- 

tipartite entanglement of the symmetry-breaking ground 




FIG. 1. (Color online) String order parameter Oz (blue) and 
staggered magnetization my (red) versus A. Their behavior 
clearly shows the existence of two different phases of the sys- 
tem: when A < 1 there is topological order protected by the 
symmetry Z2 x Z2 and a non- vanishing string order parame- 
ter, while when A > 1 there is antiferromagnetic order and a 
non- vanishing staggered magnetization. 
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FIG. 2. (Color online) Residual (multipartite) entanglement 
T in the symmetry breaking ground state versus A. From 
its behavior we can identify the two phases of the system: 
a "cluster phase" (useful for quantum computation), where 
multipartite entanglement is constantly equal to its maxi- 
mum value, and an ant i- ferromagnetic phase, characterized 
by a sharp decrease in r and the presence of a non- vanishing 
staggered magnetization along the ^/-axis. 



state is therefore able to detect the phase transition. 

The block entropy is studied in Sec. VII (see Fig. 
and yields c = 3/2 as central charge. We demonstrate 



that such a result can be traced back to the periodicity 
of the free energy, implying in turn that the system is 
energetically equivalent to three uncoupled Ising chains. 
Therefore, the criticality of CIM is characterized by an 
emergent Eg x Es x Eg symmetry [19^. 
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TABLE 1. Overview of the properties of the ground state as A varies 



Cluster state @A = 


QPT @A = 1 


Ising state @A = 00 


symmetry protected topological order 
by Z2 X Zi symmetry 


central charge c — ?>/2 
z/=z=l, /3 3/8, a = 


Z2 symmetry 



Summarizing: for A = and open boundary condi- 
tions, the system has a Z2 x Z2 symmetry that protects 
the topological order [16]; for A < 1 there is a "clus- 
ter phase" with a non-local hidden order, detected by a 
nonvanishing string-order parameter O^; for A > 1 the 
system is antiferromagnetic with long-range order. The 
antiferromagnetic and the cluster phases are separated 
by a continuous quantum phase transition. 

The critical indices of the quantum phase transition 
are summarized in Table |l] (where z and a are re- 
lated to the correlation length, the dynamical correlation 
functions and the specific heat, respectively). Finally, 
in Sec. VIII we compare our findings with other known 
models with three-spin interactions and draw some con- 
clusions. 



and Afc = Wl + A^ - 2Acos (^). Finally, we obtain 



k=l ^ ^ 



(6) 



whose ground state is defined by 7/. = 0, V/c. 

B. Free Energy 

From Eq. (|6| we can compute the partition function of 
the system: 

Z(/3,A)=Tre-^^^^^^W^^-i/^) 
= n2cosh(/3A,), 



III. THE EXACT SOLUTION 

A. Diagonalization 

Despite the presence of a three-spin interaction, the 
Hamiltonian ([T]) can be diagonalized and describes free 
fermions. Introducing the Jordan- Wigner transforma- 
tions 

'^. = fn<V7' 4=fn<V;' (2) 

\m=l / \m=l / 

where erf = (aj ± \(j^-)/2, Eq. (jlj) takes the form 

N 



= - + Q+i) 

1=1 

N 

+A^(c|+q)(c|+i-q+i), (3) 



Z=l 



apart from a border term, which is negligible in the 
thermodynamic limit. Furthermore, we apply a Fourier 
transformation bk = ZljLi ^"^^^Cj, with k = 
1, . . . , A/", followed by a Bogoliubov transformation, bk = 
Ukjk + i^/ctI/c^ with 



^ /1 , 



.±,iMh)f^^. (4) 



cos 



/47rfc\ 



Acos(^— j, 
+ Asm^— j, 



(5) 



where P = l/Zc^T with Ub being the Boltzmann constant 
and T the temperature of the system. In the thermody- 
namic limit {N 00), the free energy density reads 



/(A A) 



dplog[2cosh(/?A(p))], (8) 



where 



A(p) = a/1 + A2-2Acos(3p). 



(9) 



Notice now that, by using the periodicity of the dis- 
persion relation A(p), the free energy ([8| of the CIM ([T]) 
can be rewritten as 



1 



dp log [2 cosh {/3A^''''^{p))] , (10) 



where A^^^^^{p) = y^l + A^ — 2A cosp is the dispersion 
relation of the quantum Ising chain in a transverse field: 

N 

i^ising(A) = - ^ [ajaj+i + Xa]] . (11) 

Therefore, the system has exactly the same free energy 
of a suitable Ising model where the cluster interaction 
and the antiferromagnetic exchange along y turn into an 
effective antiferromagnetic exchange along x and an ex- 
ternal field, respectively. 

According to the general theory of continuous phase 
transitions, the second derivative of the free energy den- 
sity with respect to A is divergent: 



dimx) 



7rA2(l + A) 
- (1 + A2) K 



{\ + lf£ 



4A 



(1 + A)^ 



4A 
(T+A)^ 



cost X log|A — 1|, 



(12) 
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where £ and /C are the complete ehiptic integral of the 
second and first kind, respectively: 



S{x)-- 

lC{x)= r — = 
Jo VI 



de 

\/\ — sin^ Q 



(13) 



and the divergence of Eq. (12) in the limit A ^ 1 is a 
consequence of the singular behavior of /C(x) at x = 1. 
Such a singular behavior is ultimately due to the vanish- 
ing of the energy gap between the ground and the first 
excited state at the critical mode p = and A = 1, with 
critical indices z = = 1. 



C. Duality 




[ ♦ 




{ •I 






[as J 



The above-discussed link with the Ising model indi- 
cates that the system must enjoy a nontrivial duality. 
Indeed, consider the duality transformation 

A.,^ = cTjaJ+„ (14) 

/c=l 

(with the convention o^j^-^ — where /i^ are the Pauli 
matrices. Since the inverse of the second transforma- 
tion isa| = (with /ig = 1), one obviously gets 
from ( [lT| 

Utlt^{\) = Ai7i,ing(A-^), (15) 

apart from a boundary term. But one also gets that the 
CIM Hamiltonian ([T]) behaves in the same way: 

i^^^^^A) = Ai^(A-^). (16) 

Thus, the Ising interaction is mapped into the cluster one 
and viceversa, when the thermodynamic limit is consid- 
ered. We notice that the QPT A = 1 corresponds to a 
self-dual point [lOl [H [IS] . 

D. Majorana fermions 

An additional piece of information is unveiled when one 
expresses the CIM Hamiltonian ([T]) in terms of Majorana 
fermions (see Fig.|3|. In terms of the fermionic operators 
introduced by the Jordan- Wigner transformations ([2| , we 
have 

d2j-l = c] + Cj, d2j = i{cj - c]), (17) 
(j = l,...,iV), with 

dl = dk, {dk, di} = 25ku (18) 
(fc, / = 1, . . . , 2N). The CIM Hamiltonian becomes 

N N 

H{\) = i^d2i-2d2i^i +iA^a2^-ia2^+2, (19) 
1=1 1^1 



FIG. 3. (Color online) Each fermion of the chain (top) can 
be represented by a Majorana pair (center). The Majo- 
rana fermions can be arranged according to the interaction 
described by the Hamiltonian |T9| (bottom). The (black) 
dashed lines represent Ising interactions. The (green) con- 
tinuous lines represent the cluster interactions. The boxed 
Majorana fermions are responsible for the ground state de- 
generacy of the cluster state. 



where periodic boundary conditions have been consid- 
ered. Notice that in an open chain one gets 

N-l N-1 
H{X) = i d2l-2Cl2l+l + CL2l-ld2l+2, (20) 

1=2 1=1 

so that, by turning off the Ising interaction, four free Ma- 
jorana fermions emerge ai, as, aAr_2, ^at (Fig. [3|. This 
effect marks the non trivial Z2 x Z2 of the ground-state 
degeneracy of the system in the cluster limit A = 0. by 
tuning off the Ising interaction in an open chain four 
free Majorana fermions emerge ai, as, aAr-2, ^at (Fig.[3|. 
This effect marks the non trivial Z2 x Z2 of the ground- 
state degeneracy of the system in the cluster limit A = 0. 

Henceforth we shall focus on periodic boundary condi- 
tions. 



IV. CORRELATION FUNCTIONS 

In this section we detail the calculation of the spin 
correlation functions at temperature T, defined as: 

Rfi{T) = {cTfcTr)T with a = x,y,z, (21) 

with {■)t = ^Tr (• e~^/^^^) denoting expectation values 
in the canonical ensemble. The method we shall employ 
is a straightforward application of the techniques adopted 
in Kefs. [22, .23^. 
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FIG. 4. (Color online) Absolute value of the two-point correlation functions at fixed temperature vs A: (a) i?^(0) for r — 
3, 6, 9, 12; (b) Rr{l/2kB) for r = 3, 6, 9, 12. A = 1 is the critical point in the thermodynamic limit. 




A 



FIG. 5. (Color online) Absolute value of the two-points correlation functions at fixed temperature vs A: (a) K!f.{0) for r — 
1, 2, 6, 14; (b) Rr{3/2kB) for r = 1, 2, 6, 14. A = 1 is the critical point in the thermodynamic limit. 



Let us start considering the correlations Rji{T). We 
find 



(^i-c]) n (l-2c«4) (4+q)\ 

j<m<l I 



j<m<l 

— • • • Ai-iBi-iAi)t, 



where 



T 

(22) 



(23) 



Using Wick's theorem we can evaluate the vacuum ex- 
pectation value of a product of anticommuting operators 
Aj's and 5j's, in term of contractions of pairs thereof. 
The needed contractions are 



{BjBi)t 
{BjAi)t 



-5ji 

Dji{T) = DU-l,T) = D{r,T) 



(24) 



where the function at a fixed temperature depends 
only on the relative distance r = j — I between spins, 
as a consequence of the translational invariance of the 
system. In the thermodynamic limit the explicit form of 
this function is 



D{r,T) 



1 r 

Jo 



dp 



tanh (/3A(p)) 



X {cos [(r + 2)p] - A cos [(r - l)p]} , (25) 



where A{p) has been defined in Eq. ([9|. 

As for the standard Ising model the contractions 
{AjAj)T and {BjBj)T vanish at equilibrium. Therefore 
the pfaffians can be reduced to determinants. The result- 
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ing spin correlators, Rji = = = -R^^, read 



D{-l,T) D{-2,T) 
D{0,T) D{-1,T) 

D{r-2,T) D{r-3,T) 



D{-r,T) 
D{-r + l,T) 

D{-1,T) 



(26) 

Similarly, for the other two correlation functions we have 



RUT) = 



D{1,T) D{0,T) 
D{2,T) D{1,T) 

D{r,T) D(r-1,T) 



D(-r + 2,T) 
D{-r + 3,T) 

D{1,T) 



and 



R^T) = D{0,Tf - D{r,T)D{-r,T). 



(27) 



(28) 



Finally, we observe that the magnetization along the z- 
axis niz is given by 



m,(T) = {a|)T = (1 - 2cjc])T = {AjBj)T = -D{0,T), 

(29) 

where translational invariance has been invoked. We now 
show that the system under investigation is endowed with 
a very specific correlation pattern. To this end we notice 
that (I25I) has the form 



I)(r,T) = /(r + 2)-A/(r-l), 



(30) 



where 
I{n) 



1 r 

Jo 



dp cos(np) 

7r/3 



tanh {(3A{p)) 



dp 



1 

JO 

tanh {l3A{p)) 



cos(np) 



Hp) 

/27rn\ , / 
f 2 cos ( j cos[np) 



(31) 



X 



Hp) 



(32) 



The square bracket, and therefore /(n), is non- vanishing 
only if n is a multiple of 3, and, therefore 



r = 3m + 1, m e Z. 



(33) 



This result implies that all correlation functions in 
Eq. (28) and the magnetization along the 2:-axis (29) are 



identically zero. On the other hand, the correlation along 
the X-axis R^(T) is non- vanishing only when r is a mul- 
tiple of 3, and its absolute value at a fixed temperature 
decreases wh en the rela tive d istance increases, as can be 
seen in Figs. |4(a)| and |4(b)[ Clearly, the x-correlation 
is more sensitive to criticality for smaller r. The cor- 
relation along the ?/-axis i?^(T) is always non- vanishing 
and negative (positive) when r is odd (even), as expected 
because of the anti-ferromagnetic nature of the Ising in- 
teraction in ([1]). See Figs. 5(a) and |5(b) Its absolute 
value is small when A < 1 and tends to its maximum 
value 1 when A ^ 00. Moreover, at zero temperature 



[Fig. |5(a)| , the transition around the critical value A = 1 
becomes more abrupt for larger relative distance r, while 
at a nonvanishing temperature [= 3/2/cb in Fig. 5(b) 
the transition is softer. Finally, the curves of R^{T) be- 
come less and less peaked when distance is kept fixed and 
temperature is increased. See Fig. |6(a)[ Analogously, at 
a fixed distance, the transition to a non-vanishing value 
of R y{T) is sharper for lower temperatures, as shown in 
Fig.[6(b7 



The global picture that emerges from the study of the 
correlation functions is the following. Around critical- 
ity, Rf{T) is more peaked at small temperatures and 
small distances r; Ry{T) undergoes a sharper transition 
from zero to unity at small temperatures and large dis- 
tances r. Notice the opposite behavior of R^ and R^ for 
small/large values of r, a consequence of the order in the 
y direction. 

Using these results, in the following subsections we will 
evaluate the staggered magnetization and the string pa- 
rameter of our model. We will see that these two quan- 
tities are able to capture the nature of the Ising and the 
cluster phase, respectively, sharply marking the phase 
transition at A = 1. 



A. Staggered magnetization 

Because of the Z2 symmetry of the Hamiltonian, we 
cannot compute the staggered magnetization by directly 
applying the definition ruy = (( — l)-^crJ)T, as we would 
always obtain a vanishing result. Hence we use the for- 
mula 



lim {-iyRy{o) 



ml. 



(34) 



We notice that ( |27[ ), defining i?^, is a Toeplitz matrix, 

There- 



namely a matrix A with elements aij 
fore we can apply Szego's theorem enabling us to fix the 
asymptotics r ^ 00 of the correlation functions. Un- 
der the assumption that a^s can be considered as the 
coefficients of a Fourier series of suitable function f{p): 
am = ^ Jl^ e~^^^/(p) dp, we remind that the theorem 
states that 



^ . det I 
lim 



exp I 



(35) 



where log/i ^ J^^ log f{p) dp, and gn 

^ J^^e~^^P log f{p) dp. To apply the Szego theo- 



rem to our problem, we write 



Dir,0) = -- 



We consider the integral 



/A - e-^ip 
A-e3ip "^P- 



(36) 



D{r +1,0) = C{r) = — 



"cip) dp, (37) 
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(a) (b) 

FIG. 6. (Color online) Absolute value of the two-points correlation functions at a fixed distance vs A: (a) — (T) for T = 
0,l/2/cs,l/A:s,2/A;s,5/A;s; (b) R\^{T) for T — 0,l/2A;s,l/A^s,2/A;s,3/A;s. A = 1 is the critical point in the thermodynamic 
limit. 



with c{p) 



and 



\ r I r A - e"^'^ 

log/i = - + ^ log ^ _ ^3-^ dp = (38) 



where the second integral vanishes because the integrand 
is odd. This yields /i = — 1. With this definition we can 
write the spin correlation function along y as 



C(0) C(-l) 
C(l) C(0) 

C(n-l) C(n-2) 



C(-n+l) 
C(-n + 2) 

C(0). 



(39) 



In order to compute the ^n's, we start from the case 
A > 1. Defining a = 1/A, we obtain: 



When A < 1 one can proceed in an analogous way, ob- 
taining 



^ ngnQ-n = 



-oo. 



(44) 



n=l 



By plugging Eqs. (|38|, (|43|, and (|44|) into ([35| and ([27|) 
we obtain 



m^ = ±./lim(-l)^i?^(0) 



± (1 - A 



-2\3/8 



when A > 1 
when A < 1 



(45) 



Therefore the staggered magnetization along the ?/-axis 
captures the antiferromagnetic order in the Ising phase 
A > 1. The quantity rriy is displayed in Fig. [l] It plays 
the role of order parameter for A > 1, but it is nonlocal. 



LXJ / 

— sin(3/p). 



^=1 



As a consequence we have 



9n 



1 ^ PT^ 



sin(np) sin(3/p) dp 



ifnG5 = {3,6,9,...} 

otherwise 



2n 





and 



Finally, 



2n 



- 9 



4 ^ n 

neB 



loe^ 1 



(40) 



(41) 



(42) 



(43) 



B. String correlation functions 

We compute here the string order parameter 

0, = JimJ-l)^^aM(^n<T^j<-i<T^^ , (46) 

by exploiting the duality of the model discussed in 
Sec. IIIC The expression of the dual variables fi^ and 
/i| in l^q. (14) implies 



\k=l 



(47) 



We consider the dual ^/-correlation function between site 
j and site / 



^-1 
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which for j = 1 and I = N — 1 particularizes to 

z \ y X 



Eq. (23) we get a2j-i = Aj and d2j = iBj, so that 



V V XV 



n 



(49) 



In the thermodynamic hmit N ^ oo we can apply 
Szego's theorem and the string correlation can be cal- 
culated in terms of the staggered magnetization: 



hm {-ir{^y^y,_,)o 

N^oo 



y,dual ' 



(50) 



Because of the duality in Eq. ( 16 ), the quantity ^^^^ 
can be ca lculated with the procedure followed in Sec- 
tion IV A and is given by (45) after replacing A with 1/A, 
namely 

0,= 



',dual 



(1 - A^) 



3/4 



when A < 1 
when A > 1 



(51) 



This is the quantity plotted in Fig. [T] It is nonvanishing 
for A < 1 and captures the hidden order of the cluster 
phase. 



C. Finite Size Correlations Functions 

We conclude our analysis of correlations by briefly 
looking at the behavior of the correlation functions and 
magnetization along the z-axes at zero temperature, 
when the number of sites N is finite. We shall use numer- 
ical methods, starting from the Hamiltonian ([T]), consid- 
ering the boundary terms. In particular, we consider the 
cases N = 4, 6, 8, 10 and 12, with an exact diagonaliza- 
tion of the Hamiltonian and computing the expectation 
values (21) on the ground state of the system. We limit 



ourselves to the case of first neighbor spins. 

In contrast with the thermodynamic limit [see com- 
ments following Eq. (33)], the correlations along x and 
z and the magnetization along z are in general non- 
vanishing. However, they do vanish when the number of 
sites is a multiple of 3, that is when N = 6 and N = 12 
[see again comments following Eq. ([33|]. The absolute 
values of these "spurious" correlations decrease when the 
size of the system is increased. In particular, their max- 
imum vanishes like 0.97 x iV~^-^^, in accord with the 
thermodynamic limit. The behavior of the correlation 
along the x-axis and of the magnetization along the z- 
axis is shown in Fig. [t] (the plots of the correlation along 
z being analogous). This finite-size analysis corroborates 
the findings of the present section. 



D. Correlations between Majorana operators 

We study here the ground-state correlations between 
Majorana operators, by employing the Majorana repre- 
sentation of the system, introduced in Eq. (17). From 



{ct2j-ia2i-i)o = {AjAi)o = Sji 
{d2j-ia2i)o = i{AjBi)o = -iD{r,0) 
{d2jd2i-i)o = i{BjAi)o = iL>(-r,0) 
{d2jd2i)o = —{BjBi)o = Sji 



(52) 



where r = j — I and D{r) is defined in Eq. (25). Equa- 
tions (52) can be summarized as: 



{djai)o = Sji + i{rN) ji, j,l = l,.. .,2N, 



where 



N 



no 

Hi 



n_i 
no 



n_Ar+i 
n_Ar+2 



njv_i nAr_2 ... no 

D(i,0)" 
-D{-j,0) 



(53) 



(54) 



(55) 



with D(r, 0) defined in Eq. (25). These expressions 
will be useful in the following, in particular in Sees. [V| 
andlVHl 



V. REDUCED DENSITY OPERATORS 

We provide here explicit expressions of the reduced 
density matrix for a single spin, two spins and a block of 
L spins. Let us briefly derive the reduced spin density 
matrices by suitable partial tracing of the state p of the 
spin system: 

{aj^^^ ...aj^\n^){Q^\aj^^^ ...aj^). (56) 

where \^w) denote the states of the chain arising with 
probability and a G {t^i}- Here (ji,...,jAr) is a 
given permutation of (1, . . . , TV). 

First we consider the reduced density matrix of a block 
of contiguous spins; the single and two spin cases will 
be obtained as particular examples. By translational in- 
variance we can consider, without loss of generality, the 
block of the first L spins. Then we can expand the matrix 
P{L} = P12...L as 

PW=2"'^ E (57) 

Q;i,...,Q;i,G{0,a:,2y,2} 



where 



(58) 



can be obtained through the expressions of the Majo- 



rana operators outlined in Section IV D In fact, we know 
that the system is invariant under parity transformation. 
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hence Pai...aL — whenever the sum of the a's equal to 
X and of the a's equal to y is odd. Therefore the non- 
vanishing coefficients of the expansion (57) correspond 



to the expectation values of a product of Pauli matrices 
with an even total number of cr^^'s and a^y's. Products of 
this kind are mapped into products of an even number 
of Majorana operators a^, with j = 1, . . . , 2L. We can 
then can use the Wick theorem to express such products 
in terms of the correlations {ajdi) which are contained in 
the matrix F^r of Eq. (53). 



Let us now consider the state of a single spin placed at 
site j of the chain. In general, the spin will be in a mixed 
state. The Hilbert space is and in the basis {cr^ = 
Ij ^(j^ ^aj} the single-spin reduced density matrix is 
given by 



Pj 



ie{0,x,y,z} 



(59) 



where Pa = Tr {cr^ pj^ = {cr^)- The system considered in 
this article is traslationally invariant, and therefore the 
single-spin state, that we denote pi, is translationally 
invariant as well. 

Finally, for the case of two spins placed at sites j and / 
of the chain, the matrix pji acts on the Hilbert space C^(8) 

of which the set |cr^ |, with a, /3 G {0,x,?/,2:}, 

is an orthonormal basis. Therefore we can make a formal 
expansion of the reduced density operator on such a basis 



Pjl = -A 



a,(3E{0,x,y,z} 



(60) 



VI. SPIN ENTANGLEMENT IN THE GROUND 
STATE 

In this section we analyze the entanglement between 
two spins and between a single spin and the rest of the 
system. As measures of such quantity we will make use of 
concurrence [24] and residual entanglement [25 respec- 
tively. The concurrence is defined as 



C(pi,.+i) = C(r) 



(61) 



where 71 > 72 > 73 > 74 > are the eigenval- 
ues of the matrix IZr = Pi,r+iPi,r+i, with pi^r-\-i = 
{gy (g) gy) (^ ^ ^ ^^)? where pi,r+i is the reduced 



density matrix (60) and the complex conjugation is taken 



in the computational basis. Concurrence varies from 
C = for a separable state to C = 1 for a maximally 
entangled state and is a convex functional: 



The residual entanglement (tangle) is defined as 



N-l 



T = 4det/9i - C^{r) 



(62) 



(63) 



r=l 



where pi is the 1-spin reduced density matrix (59). The 
tangle is a global measure of the multipartite entangle- 
ment encoded in a (pure) state. 



where Pa (3 = Tr ^a^af pj/^ = {g^g^). Equation (60) 

makes it manifest that the state of the couple of spms 
j and / can be written in terms of the spin correlation 



functions computed in Sec. IV 



A. Entanglement of the "thermal" ground state 

Let us consider the concurrence between couples of 
spins when the system is in the "thermal" (completely 
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mixed) ground state po arising as the limit (3 
thermal canonical state p = (l/Z) e~^^: 



Po 



(|(^l)(l]l| + |(^2)(^2|), 



oo of the 



(64) 



where |^^i,2) are the two degenerate ground states (we 
remind that we are considering periodic boundary con- 
ditions). By construction, po has the same symmetries 
as the Hamiltonian and, therefore, cannot spontaneously 
break any symmetry. 

A straightforward computation yields the following 
eigenvalues of the matrix IZr 



71 = 






+Rmf 


72 = 




+ J??;(o) 


- Rmf 


73 = 




- i?^(0) 


+Rmf 


74 = 




- i?^(0) 


- Rmf 



(65) 



We notice that, when r is not a multiple of 3, the correla- 
tion vanishes, so that the eigenvalues become pairwise 
equal: 



71 = 73 = ^ (1 + Rmf 

72 = 74 = ^ (1 - Rmf 



(66) 



This implies that the concurrence C{r) vanishes when r 
is not a multiple of 3. On the other hand, when r is a 
multiple of 3 the four eigenvalues are in general distinct 
from each other but the difference of their square roots 
remains negative, making the concurrence vanish also in 
this case. Figure [S] shows the behavior of the eigenvalues 
for r = 3 and r = 6. 



It is instructive to see how the concurrence vanishes in 
systems of small size. We numerically analyzed this prob- 
lem and our main results are shown in Fig. [9] In the ther- 
mal case, when = 4, the concurrence is non- vanishing, 
both for first and second neighbors; when N = 6 the 
concurrence always vanishes, while for N = S only the 
first-neighbor concurrence is non- vanishing, with a value 
that is much smaller than that for = 4; finally for 
A" = 10 and A" = 12 the concurrence always vanishes, as 
in the thermodynamic limit. The particular behavior for 
N = 6 can be viewed as a consequence of the fact that 
the correlations have the same value as in the thermody- 
namic limit, as seen in Sec. |IV C[ For A" = 4 and 8 we 
have finite-size effects similar to those seen in the case 
of the correlation functions. The case A' = 10 is partic- 
ularly interesting: even though the correlation functions 
behave differently from their thermodynamic limit (see 
Fig. [t]), concurrence identically vanishes. The thermo- 
dynamic limit appears therefore to be reached earlier in 
terms of concurrence than in terms of correlation func- 
tions. 



In conclusion, for the "thermal" ground state (64), 

C(r) = Vr. (T = 0) (67) 

We conclude that in the "thermal" ground state there is 
no entanglement between any pairs of spins, at any dis- 
tance from each other, for any value of A. The concur- 
rence of the "thermal" ground state is therefore unable 
to signal the phase transition at A = 1. 

Incidentally, we observe that the same analysis can be 
applied to the canonical thermal states at finite temper- 



atures, by replacing R^{0) with R^{T) in Eqs. (|65f . The 
eigenvalues 7^'s have the same properties reported above 
and therefore, as was to be expected, the result (67) is 
valid at any temperature 

C{r) = Vr. (T ^ 0) (68) 

Finally, let us compute the residual (multipartite) en- 
tanglement r, defined in Eq. (63), when the global system 
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FIG. 10. (Color online) First-neighbor concurrence vs A, for 
FIG. 9. (Color online) Non- vanishing concurrences vs A in different values of the intensity h of the symmetry-breaking 
the "thermal" ground state of finite-size systems. staggered magnetic field and iV = 8. The points for /i = are 

taken from Fig. [9] (A^ = 8). 



is in the "thermal" ground state (64). One obtains 



r = 1, VA. 



(69) 



Eqs. (67) and (69) enable us to conclude that neither the 



two-spin nor the multipartite entanglement of the "ther- 
mal" ground state are able to detect the quantum phase 
transition. A natural question arises: does the entangle- 
ment of the symmetry-breaking ground state detect the 
phase transition? This problem will be tackled below. 



B. Entanglement of the ground state with broken 
symmetry 

We now discuss the effect of spontaneous symmetry 
breaking on the ground-state entanglement. Symmetry 
breaking is achieved by adding to the Hamiltonian a little 
staggered magnetic field along the y-axis h^j{—iya^^ 
that breaks the invariance of the Hamiltonian under par- 
ity transformation. We start with two-spin entangle- 
ment (concurrence), by performing a numerical analysis 
on small systems, like in Sec. |IV C 



We start by observing that the convexity of the concur- 
rence, Eq. (62), when applied to the "thermal" ground 



state (64) yields 



= C (po) < 2 [C (Pi) + C (P2)] = C (pi) , (70) 



due to (67) and the symmetry between pi = 
and p2 =1^2) (^2 1- Therefore, in principle, the concur- 
rence of the symmetry-breaking ground state could be 
nonvanishing (and possibly detect the phase transition). 
We now show that this does not happen. 



Figure 10 displays the first-neighbor concurrence for 
the symmetry-breaking ground state of a finite chain of 
N = S spins and a varying strength of the staggered 
magnetic field. The plot indicates that concurrence is 



negligibly affected by the symmetry breaking mechanism: 
indeed, for ^ the curve becomes identical to that of 
Fig. [9] at = 8 (approximating it from below). The 
concurrence of second and third neighbors, as well as 
for longer chains {N > 10) yields identical results. This 
is an evidence that there is no bipartite entanglement 
between couples of spins of the chain, even when one 
of the degenerate ground states is chosen, breaking the 
parity symmetry. Therefore we conclude that two-spin 
entanglement is unable to signal the phase transition. 

We finally analyze the effect of symmetry breaking on 
the residual entanglement of the ground state. To this 
end, we compute the state of spin j when one of the 
degenerate ground states is chosen. 

The expansion of Sec. |VI A] is still valid, but the values 
of the coefficients Pa [that define the state, see Eq. (59)] 
are modified. The quantities that are symmetric under 
parity transformation are unaltered, and therefore the 
magnetization along the z-axis is still vanishing. On 
the other hand, the values of (aj) and (aj) can be fi- 
nite. With a calculation very similar to that presented 



in Sec. IV A one easily shows that 



{a^)l = lim Rm 



(71) 



vanishes for any value of A. Therefore the reduced density 
matrix of spin j is 



1 



(72) 



On the other hand the concurrence vanishes and therefore 
the residual entanglement is 



4detpj 



1 



1 - (1 - A- 



.3/4 



if A < 1 

if A > r 

(73) 
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Block Length (L) 



FIG. 11. (Color online) Block entropy Sl (crosses) and its 
logarithmic least square fit (full line) versus the size of the 
block L. 



This is shown in Fig. [2j This result indicates that the 
multipartite entanglement of the ground state with bro- 
ken symmetry saturates in the cluster phase (A < 1), 
whereas it decreases in the antiferromagnetic phase (A > 
1). The quantum phase transition is detected by a sin- 
gularity of its first derivative: 



dxr{X) 



-cost X (A — 1) 



-1/4 



(74) 



VII. BLOCK ENTROPY 



In this section we compute the entanglement of a block 
of contiguous spins. As a measure of entanglement, we 
take the Von Neumann entropy: 



Sl = Tr {p^Ly log^ P{l}) , 



(75) 



where p^^y is the density matrix of a block of length L 
defined in Eq. (57) and does not depend on the position 



of the block because of the translational invariance of the 
system. 

To calculate 5'l we will follow the procedure originally 
developed in Ref. [27 . It results 



5l = E^'''^'^ 



(76) 



where 



H{x) = -X log2 x-{l-x) log2(l - x) (77) 



is the Shannon entropy, and ii^^j, (j = 1, . . . , L) are the 
(purely imaginary) eigenvalues of the Majorana correla- 
tion matrix Tl in Eq. (53). 

We have evaluated Sl numerically with blocks of 
length ranging from 2 to 200 spins, at the critical point 



A = 1. The resulting behavior, displayed in Fig. [TT] is 



6'l - log2 L^a, 



(78) 



where the least square fit yields for the constant multiply- 
ing the logarithmic term the value 0.506 with a standard 
error 0.001, while a ^ 1.236 with an error 0.008. This 
value of a is in agreement with that one would obtain by 
self-dual arguments 



l + 7^ + (2-6/3)ln2-ln3 
21n2 



1.246, 



(79) 



where — 0.577 is the Euler-Mascheroni constant and 
/3 0.022 another constant [39 . 

The logarithmic scaling of the entropy at the critical 
point is a result obtained also in other models [28; . More- 
over the multiplicative constant of the logarithmic term 
is related to the central charge of the 1+1 dimensional 
conformal theory that describes the critical behavior of 
the chain, trough the relation [26^ 



Sl 



6 



■ log2 L, 



(80) 



where c and c are the central charges of the so-called 
holomorphic and antiholomorphic sectors of the confor- 
mal field theory. Also the central charge is directly re- 
lated to the universality class of the system, so that, if 
two quantum one-dimensional models belong to the same 
universality class, they have the same central charge. 

The result (78) shows that the central charge of our 
critical chain is c = c = 3/2. See Table |lj This value is 
different from that of the Ising chain, where cising = 1/2, 
so that our model is in a different universality class than 
the Ising model in transverse field. 

Below we exploit the periodicity of the free energy (see 
Sec. Ill) to prove that the CIM Hamiltonian can be in- 



deed recast in a sum of three decoupled Ising Hamiltoni- 
ans. From Eq. ([6|, by taking N = 3M and splitting the 
sum in three parts, we get 



H{\) 



Ising' 



where 



M 



k=l 



(81) 



(82) 



with 7^''^ = 7fe+(5-i)M and 



Ar' = w'l + A»-2Ac„.(?g:), 



(83) 



Recalling that the central charge is an extensive quantity 
[38] , we conclude that the criticality of the CIM arises as 
a simple sum of three Ising central charges. 
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VIII. CONCLUSIONS AND OUTLOOK 

We studied the statistical mechanics of the cluster- 
Ising Hamiltonian ([T]), providing important elements for 
the understanding of the phase diagram of the system. 
The specific correlation pattern that we found indicates 
that most of the physics of the system is displayed be- 
yond two-spin correlations. This is manifest in the type 
of entanglement encoded in the ground state that is of 
multipartite form. The QPT occurring at A = 1 is exotic 
because the two phases enjoy distinct symmetries that 
cannot be continuously joined; the cluster phase does not 
admit any local order parameter. 

The universality class of the Cluster-Ising model ([T]) is 
summarized in Table |lj Below we compare our findings 
with other spin models with triplet interaction in (1+1) 
dimension enjoying a Z2 x Z2 symmetry and displaying 
a continuous phase transition [29] [30] . The current un- 
derstanding of such class of models benefits of the ex- 
act solution by Baxter and Wu [ST, corroborated by the 
calculation of the spontaneous magnetization by Joyce 
[32]. Their critical properties can be extracted through 
the analysis of the quantum criticality of quantum spin 
chains with mixed two and three spin interactions in an 
external field [33 . There is general consensus that the 
universality class of such a class of quantum models at 
zero external field is identified by the criticality of the 



four-state Potts model [34[ |35] at equilibrium: v = 2/3, 
13 = 1/12, with some controversy on the dynamical crit- 
ical exponent 2; ~ 2 — 3 [36], and central charge c = 1. 
When the external fields are switched on, the degeneracy 
of the ground state results to be affected in a nontrivial 
way. Here we mention that for a transverse fields, and for 
moderate longitudinal field, the resulting models are still 
in the universality class of the four-state Potts model; for 
higher longitudinal fields, criticality is shared with the 
three-state Potts model, yielding c = 4/5 [33^,^37. The 
criticality of all these models appears to be very different 
from the Cluser-Ising model investigated in the present 
work. 

Finally, it would be interesting to study this model 
from an experimental point of view, for instance in the 
context of adiabatic quantum computation. In this case, 
one could study how the encoding of the cluster phase 
starting from the antiferromagnetic one is affected by the 
critical point at different temperatures. 
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